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Harmonic moments are integrals of integer powers of z = x + iy over a domain. Here the domain 
is an exterior of a bubble of air growing in an oil layer between two horizontal closely spaced plates. 
Harmonic moments are a natural basis for such Laplacian growth phenomena because, unlike other 
representations, these moments linearize the zero surface tension problem (Richardson, 1972), so 
that all moments except the lowest one are conserved in time. For non-zero surface tension, we 
show that the the harmonic moments decay in time rather than exhibiting the divergences of other 
representations. Our laboratory observations confirm the theoretical predictions and demonstrate 
that an interface dynamics description in terms of harmonic moments is physically realizable and 
robust. In addition, by extending the theory to include surface tension, we obtain from measure- 
ments of the time evolution of the harmonic moments a value for the surface tension that is within 
20% of the accepted value. 
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I. INTRODUCTION 



A. Laplacian growth and viscous fingering 

Non-equilibrium processes give rise to a variety of pat- 
terns with remarkable geometrical and dynamical prop- 
erties [TJ |2]. Often these processes are represented by 
the dynamics of an unstable interface between different 
phases, and the interface patterns can exhibit univer- 
sal features [3]. Examples include crack propagation [2], 
fluid- fluid interface dynamics [4], crystal formation [5], 
and biological growth [5J. 

The simplest process leading to unstable universal pat- 
terns is Laplacian growth, where the velocity of an inter- 
face is proportional to the gradient of function that is 
harmonic outside (or inside) the interface. Here we ex- 
amine the simplest example of Laplacian growth: quasi- 
two-dimensional (2D) viscous fingering in a Hele-Shaw 
cell [7J, where a viscous fluid is displaced by an invis- 
cid fluid between two horizontal closely spaced parallel 
plates. 

Figure [T] shows four viscous fingering patterns grown 
in the radial Hele-Shaw cell described in Section llVl Vis- 
cous silicone oil is removed from a buffer surrounding the 
plates and air enters between the plates through a central 
hole in the bottom plate. 

As the air bubble expands, the oil/ air interface is un- 
stable. The depth averaged velocity v(x, y) and the pres- 
sure p(x,y) in the oil are approximated by Darcy's law 
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FIG. 1: Viscous fingering patterns of air in oil contained be- 
tween two plates separated by 125 /im (image size is 14 cm x 
14 cm) (a) Bubble grown at a constant pumping rate of 0.52 
mL/min. (b), (c), (d) Bubbles grown with varying pumping 
rates, as described in Section IV. B. The time development of 
bubble (b) is shown in Fig. 3. 



where b is the spacing between the plates and /i is 
the dynamic oil viscosity. The oil is incompressible, so 
divv = 0. From (JlJ it then follows that in the oil 

V 2 p = . (2) 



'Electronic address: swinney@chaos.ph.utexas.edu 



Since the normal velocities of the interface V and of the 
fluid at the interface coincide [5], 
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where d n is the normal component of the gradient. Be- 
cause the air is nearly inviscid, the pressure in the air is 
essentially uniform, and its value can be taken as zero. 
Thus the pressure jump across the oil/air interface coin- 
cides with the oil pressure at the interface and is given 
by [J 
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where a is surface tension and k is the local curvature 
of the interface in the horizontal plane. (An updated 
boundary condition including a 7r/4 multiplicative fac- 
tor and an additional term correcting for wetting will be 
presented in Section IIIB ) 



The asymptotic pressure boundary condition in oil far 
from the interface in the radial geometry is 
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log(x 2 + y 2 ) (for x 2 + y 2 — >• oo), (5) 



where Q is the pumping rate (the rate of area growth), 
which may depend on time. Equations (JT])-(5]) complete 
the description of 2D viscous fingering, which is a proto- 
type of Laplacian growth. 

Laplacian growth, a classical free-boundary problem 
[9], is famous in both physics and mathematics: In 
physics, because the process is highly unstable, dis- 
sipative, non-equilibrium, and nonlinear, and because 
the growth produces universal patterns |H [THl EG]. In 
mathematics, because the description ([T])-(|5]), as simple 
as it looks, reveals a powerful and profound structure 
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Laplacian growth occurs in many physical systems and 
is known by names such as crystal growth, amorphous so- 
lidification [5], electrodeposition [T5J[[B], bacterial colony 
growth [B] , diffusion-limited aggregation (DLA) [T7] , mo- 
tion of a charged surface in liquid Helium [TS], and vis- 
cous fingering [19]. There are thousands of articles (the- 
oretical, experimental, and computational) devoted to 
Laplacian growth if one includes work on closely con- 
nected problems such as the Stefan problem (solidifica- 
tion) [13 EE], DLA [22], and the phase-field model [2J. 
Nevertheless, many features remain unexplained, despite 
extensive effort and full knowledge of the laws of physics 
describing the process. 



B. Harmonic Moments 

The present work concerns a powerful description of 
Laplacian growth called harmonic moments. Our work 
demonstrates that robust results for the time evolution 
of harmonic moments can be obtained from viscous fin- 
gering data. Harmonic moments are defined as 
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where k = 0,1,..., oo, z — x + iy, and the domain of 
integration D is exterior to a pattern's boundary T and 
bounded by a large circle on the outside. These are exte- 
rior moments, which are relevant for exterior Laplacian 
growth, when a viscous fluid is outside the interface, as 
in this work. Interior moments, which arc relevant for in- 
terior Laplacian growth when a viscous fluid is inside the 
interface, are defined similarly, but with positive powers 
of z under the integral, which is taken over the interior 
of the pattern. 

Harmonic moments are of fundamental importance for 
Laplacian growth because in the absence of surface ten- 
sion, there is an infinite number of conservation laws: all 
moments Mk (except for k — 0) are conserved in time, 
as was discovered by Richardson (1972) (33] for the in- 
terior Laplacian growth problem. The conservation laws 
have been extended to the exterior case [24] and are ex- 
perimentally confirmed in this work. In the harmonic 
moments basis the whole evolution of D(t) is reduced to 
the time-dependent area of a growing bubble, which is 
Mo(t). This result selects the 'harmonic basis' Mk from 
all other bases as the best basis for describing D(t) in 
Laplacian growth. 

Harmonic moments form a complete basis for repre- 
senting any 2D interface, regardless of its complicated 
shape, provided that D is analytic and singly connected 
|25l I26j . Like Fourier modes, each harmonic moment 
Mk corresponds to a particular aspect of the interface. 
The moment M is the area of the domain D divided by 
7r. Moments Mk for k > 1 are proportional to the am- 
plitude of a monochromatic wave, which modulates the 
circle with an exactly k waves along the circumference, 

z = R exp(i</>)[l + etfc exp{— ik(p)], 

where <j) is a 'stream function' parameter along the in- 
terface, ak is an amplitude, and R is the bubble radius. 
Specifically, Alk = a,kR 2 ~ h . For an approximately 71- 
fold pattern, the dominant moments are M n , Mi n , M% n , 
etc. [e.g. Fig. Ejjb) is approximately 7-fold symmetric 
and Fig. [TJc) is approximately 5-fold symmetric]. How- 
ever, care should be taken in comparing different mo- 
ments since they have different units; the moment Mk 
has units cm 2_fc . 

Viscous fingering structures become extremely com- 
plex for high growth rates, where tips repeatedly split 
and form new fingers and fjords. Hence it is remark- 
able that theory predicts that a set of purely geometrical 
quantities, Mk, will change slowly during the growth pro- 
cess and will have a well-defined limit as the surface ten- 
sion approaches zero. The unique properties of harmonic 
moments have been already used to establish connections 
between Laplacian growth and other fields of physics and 
mathematics |T3J [T3J [2] , but no previous studies have at- 
tempted to extract harmonic moments from laboratory 
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data. Our work demonstrates that robust moments Mk 
can indeed be obtained from experiments, and particu- 
larly that the moments Mf. for k > 1 are conserved in 
the limit of zero surface tension. 



C. Overview 

In the following section we present some additional 
properties of harmonic moments. Section III extends the 
theory for harmonic moments in viscous fingering to the 
physical case, where the surface tension is not zero. In 
this case, the moments Mk with positive k are no longer 
conserved, but their time derivatives are proportional to 
surface tension and hence vanish in the zero surface ten- 
sion limit. Section IV presents our experimental and data 
analysis methods. Section V presents results for har- 
monic moments determined from experiments for a wide 
range of conditions. Also we show that the results for 
the harmonic moments can be used to deduce a physical 
parameter, such as surface tension. Section VI discusses 
the significance of our results. 

II. HARMONIC MOMENTS: HISTORY AND 
APPLICATIONS 



up and down infinitely, the moments are defined as M/~ = 
J D exp(— 2irkz/L) dxdy/ir. This geometry is relevant for 
Laplacian growth in a rectangular Hclc-Shaw channel, 
which has been studied since the work by Saffman and 
Taylor g|. 

B. Connections and Applications 

In science and engineering it is often of interest to find 
the shape of an object from the indirect measurement of 
the harmonic moments Mk- The problem of domain re- 
construction from its moments has applications in many 
areas, including signal processing, probability and statis- 
tics, tomography, and the inverse potential problem in 
geophysics (magnetic and gravitational anomaly detec- 
tion) 28J. While the recovery of shapes from moments 
is often an ill-posed problem, it was recently recognized 
that the moments problem allows the complete closed 
form solution for so-called quadrature domains [215], a 
branch of mathematics created in 1970s [SOJ EH- Re- 
markably, these solutions are based on a technique of nu- 
merical linear algebra that yields numerically stable and 
fast algorithms and exposes a deep connection between 
harmonic moments and a theory of analytic approxima- 
tion [321133]. 



A. The Inverse Potential Problem 

The origin of harmonic moments dates back to Isaac 
Newton's study of the Inverse Potential Problem, which 
was subsequently investigated by Kelvin, Poincare and 
many researchers in the twentieth century. The Inverse 
Potential Problem 25 , 26J asks how to recover the shape 
of a domain occupied by a uniformly distributed mass, 
given the gravitational potential created by this mass 
outside the domain. Surprisingly, this classical math- 
ematical problem was found to be at the heart of the 
underlying mathematical structure of Laplacian growth 
[T2l IT31 H¥l |2"5] , In two dimensions the gravitational po- 
tential created by uniformly distributed mass (with unit 
density) occupying the domain D is given by 

&(x, y) = — [ log \z — z'\ dx'dy' = 
k Jd 
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if z and the origin lie outside D. The Mk, defined 
in ([6]), are multipole moments of the mass distribution 
[2"T] . In Laplacian growth the Mk are traditionally called 
harmonic moments. In ^ the moments Mk are de- 
fined for a radial geometry. For a rectangular geome- 
try, where matter occupies a horizontal semi-infinite strip 
D with a width L bounded from the left by an arbi- 
trary curve, and periodically extended (repeated) both 



III. THEORY OF HARMONIC MOMENTS 
WITH NONZERO SURFACE TENSION 

A. Derivation of dM^/dt 

In Laplacian growth with nonzero surface tension the 
harmonic moments are no longer conserved. To relate 
the dynamics of the harmonic moments to measurable 
quantities in viscous fingering, we now obtain an expres- 
sion for the time derivative of the harmonic moments in 
terms of a line integral over the air/oil interface. The 
time derivative of Mk follows from the definition |6| and 



(1) 



dMk 
dt 



dl 



We want to exchange the derivative between p and z k , 
so we add and subtract pd n z~ k from the integrand, 



= Lz- k d n p-pd n z- k +pd n z~ k )-. (8) 

Here I is the arc length along the interface T. The integral 
of the first two terms is zero by Gauss's theorem, after 
using (|3j) and the analyticity of z~ k in D, 



(z- k d n p-pd n z- k )dl = <h (z~ k Vp-pVz' k ) -ndl 
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Using Q for the pressure, we have 



dM k 
dt 
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with the interfacial curvature given by n = d6 /dl, where 
9 is the angle of a tangent line to the interface. Using 
the identity d n z~ k — idiz~ k = —ikz~( k+1 ^dz/dl, where 
di denotes a tangential derivative, we have 



although in these experiments there is small motion at 
the base of the fjords as a consequence of relaxation due 
to surface tension; the base of the fjord becomes more 
bulbous the longer surface tension has acted on it. 

After substituting this new expression for the pressure 
jump into (10 1, we obtain 
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Testing the theory and determination of surface 
tension 
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where the last expression was obtained using dz/dl 
exp(i#). Integrating by parts, we obtain 



(12) 



Note that Richardson's result for the conservation of mo- 
ments is recovered in the limit that the surface tension a 
vanishes. Surface tension is now a regular rather than a 
singular perturbation. 



B. Wetting and scaling correction 

The pressure jump at the interface Q must be modi- 
fied because as the air bubble advances between the glass 
plates, it leaves behind a wetting film on each plate. The 
thickness of this film increases with the interface's ve- 
locity [34 . This film can be seen as interference fringes 
on the images in Fig. [I] because the interface is slowing 
down as the bubble expands. In addition, a factor 7r/4 is 
needed in Q. The pressure jump including the wetting 
correction and a factor 7r/4 [not in Q] was calculated by 
Park and Homsy to be (34] , 
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where V is the local normal velocity of the interface. The 
thickness of the film predicted by Park and Homsy was 
experimentally confirmed by Tabeling and Libchaber [35 
for 6 x KT 4 <^<3x 10~ 3 . 

(7 

The interface velocity V is greatest at the finger tips 
and is much smaller at the sides of the fingers. Conven- 
tionally the base of the fjord is called a stagnation point, 



In Section V we compute the moments Mf.(t) for grow- 
ing viscous fingering patterns directly from the pattern 
geometry using the definition (6j). The time evolution of 
the moments can be compared to the result in the zero 
surface tension limit, where the moments are conserved. 

We will also compare different moments by computing 
the normalized amplitudes, 



a k = \M k \/\M \ 



(15) 



Further, we can directly test the theory using ( 14) and 



measuring the time derivative of the moments and inde- 
pendently determining: 

i. the fluid surface tension a, 

ii. the fluid viscosity /z, 

iii. the thickness b of the gap between the plates, 

iv. the interface velocity V in the second integral in 
p"4| (this can be obtained from analysis of images 
of the air bubble as a function of time, described in 
Section 
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v. the two integrals in (14) 



Alternatively, we can use the theory to determine a 
fluid property if the other four quantities arc indepen- 
dently known. This is the way in which we determine 
the surface tension in Section V. However, since numeri- 
cal differentiation of the data for Mk(t) is difficult to do 
accurately, we take one more step to obtain a working 
equation with improved signal to noise: equation (14 1 is 
integrated over a time interval to obtain 
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With only the definition of moments (|6]) , the dynamics 
cannot be related to known physical quantities. However, 
with (16), the dynamics can be quantitatively connected 
to known physical quantities. Equation 16 is a cubic 
equation of the form aa + ba 1 ^ + c = 0, where a, b, and c 
are complex numbers. We solve this equation numerically 
in Section V.B to deduce the surface tension and to test 
the predicted decay rates of the moments. 



IV. EXPERIMENT 
A. Apparatus 

An oil layer is contained in a Hclc-Shaw cell consist- 
ing of two horizontal, closely spaced glass plates with a 
hole through the center of the bottom plate. When oil is 
pumped out of a buffer that surrounds the oil layer, air 
(at nearly atmospheric pressure) enters the layer through 
the hole in the bottom plate and forms a bubble in the 
center of the oil layer. 

The optically polished glass plates each have diameter 
28.8 cm and thickness 6.0 cm; each plate is flat to 0.2 
/itm, as described in [TU1 US] • The gap between the plates 
was either 125 ± 5 /im or 384 ± 6 /zm; most of the gap 
uncertainty arises from using a micrometer to measure 
the thickness of the metal shims that set the size of the 
gap. Interferometric measurements using a sodium lamp 
showed that the gap thickness was uniform to 0.3 ^m for 
the 125 fim. gap and 1.6 /mi for the 384 ^m gap. 

The oil was Dow Corning 200 silicone oil at 24 °C. 
We measured the viscosity /x = 49.9 ± 0.3 mPa s with a 
Paar Physica MCR300 rheometer. The surface tension 
^reference = 21.1 ± 0.1 mN/m was measured by the Wil- 
hclmy Plate method using a Kruss Kll tensiometer. The 
density was measured to be p — 0.9585 ± 0.0005 g/cm 3 . 

Bubble patterns were imaged from above with a CCD 
camera (1300 x 1030 pixels). The frame rate ranged from 
1/6 to 2 frames/s, depending on the growth rate of a 
bubble. 



B. Growth of a bubble 

Experiments were initiated by obtaining a nearly cir- 
cular air bubble, grown by slowly withdrawing oil from 
the buffer surrounding the gap between the plates. Af- 
ter an initial nearly circular bubble had been grown to 
a radius of at least 2 cm (to give good spatial resolution 
in the images), multi-fingered bubbles like those in Fig. 
1 were grown by using a precision computer-controlled 
syringe pump to remove oil from the annular buffer at a 




0.10 

pumping rate q (mL/min) 



FIG. 2: In order to grow an approximately n-fold symmet- 
ric bubble, the pumping rate is adjusted to follow (17 1, as 
illustrated here by curves for different n. The 7- fold bubble 
in Fig. 1(b) was grown in this way, as indicated by the bold 
line. We measure the area of a growing bubble in real time 
from the video images, and compute from that area a radius 
for a circle with that area. Then, given a target n-fold sym- 
metry, the pumping rate is adjusted (in small steps) to stay 
on the desired n-fold curve of the graph. (These curves apply 
to a cell with gap thickness of 125 (im.) 



specified rate. Growth of a full-sized bubble (15-20 cm 
diameter) took from 30 to 1600 s (typically 300 s). 

After obtaining a nearly circular bubble, a multi- 
fingered bubble can be obtained by pumping oil out of 
the buffer. Usually fluid is removed at a fixed rate, as in 
our laboratory's previous viscous fingering experiments 
|10[ 119] . However, with a fixed pumping rate, the n-fold 
mode that is the fastest growing changes with time, and 
the resultant bubble has many azimuthal modes with a 
substantial amplitude [e.g. Fig.[ljd)]. 

Our procedure for obtaining the harmonic moments 
and deducing the surface tension is applicable for bub- 
bles grown with any pumping rate, q{t) = bQ(t). How- 
ever, results for the harmonic moments are more robust 
if a bubble retains an approximate n-fold symmetry as it 
grows. To achieve this, the pumping rate must decrease 
in a way that maintains an approximate n-fold symmetry. 
A linear stability analysis by Bataille [36j [37] provides an 
expression for the pumping rate needed to maintain the 
same n-fold perturbation of a circle as the fastest growing 
[Ml IS?] mode, 



Qn = 



6pR 



(3n 2 - 1) 



(17) 



where R is the radius of the bubble. The bubble radius R 
given by the Bataille formula is plotted as a function of 
the volumetric pumping rate q = bQ for different n-fold 
modes in Fig. [2] Conventional pumping with constant Q 
would appear as a vertical trajectory in Fig. [2] For fu- 
ture experiments one should consider using Park-Homsy 
pressure jump, p = —amr/i, instead of using Bataillc's 
assumption that p — —an. 

The particular multi-fingered pattern that develops de- 
pends not only on pumping rate but also on initial pertur- 
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bations and non-uniformities of the glass plates. In the 
early stages of the pattern growth, when a bubble begins 
to deviate from a circle, small unavoidable departures 
from the circular shape have a greater influence on the 
development of the interface than the accuracy of pump- 
ing or non-uniformities of the glass plates. In practice, 
we found the smallest m achievable was 5; for smaller 
target m the growth was too slow to overcome initial 
perturbations on the bubble's surface before the bubble 
reached the edge of the cell. We found that although 
the application of the Bataillc formula for small pertur- 
bations of a circular bubble is not justifiable for bubbles 
that have grown large fingers, in practice the Bataillc's 
expression remains helpful in maintaining a symmetrical 
target pattern. 

The present work concerns bubbles that are growing 
throughout an experiment. If the pumping were stopped, 
the surface tension would begin to smooth the interfacial 
regions of high curvature, and an interface could even re- 
verse direction, absorbing the wetting film it left behind. 

In most viscous fingering experiments where oil is with- 
drawn at a constant rate, the process has a natural sep- 
aration of time scales: tq — A/Q, determined by pump- 
ing rate, is usually much less than the "capillary" time, 
t a = 48fiR 3 /(crb 2 n(n + 1)) (for not very high n), which 
corresponds to smoothing of a pattern by surface ten- 
sion. The difference between time scales (t a >> tq) 
makes possible the rich interfacial patterns. However, 
in our case there is only a single time scale because the 
pumping Q is adjusted continuously to maintain approx- 
imate n-fold symmetry. Equating the scales £q ~ t a , one 
recovers (within a multiplicative factor that depends on 
n) the Bataillc formula. Note that for constant pumping 
the radius of a bubble is given by R ~ y/i, while in our 
case R ~ \ft (because Q ~ 1/R and R 2 ~ Qt). 



C. Image analysis 

The oil/air interface in each image was first obtained 
by subtracting the background image, thresholding, and 
using an edge detection algorithm. This located the in- 
terface to within a pixel. The resolution of the interface 
was typically 50 pixels/cm, so that one pixel's length was 
about 0.2 mm. 

The position of each point on an interface was then 
obtained with sub-pixel accuracy by interpolating the lo- 
cation along a line perpendicular to the rough interface 
that was half of the intensity difference between the in- 
side of the bubble (more intense) and the outside of the 
bubble (less intense) . The algorithm typically found a po- 
sition of half-intensity to 0.1 pixel (about 20 fim, which 
is small compared to the 125 /zm or 384 /im plate sep- 
aration). This procedure yielded a sufficiently smooth 
interface so that smoothing was not needed. 



D. Determing surface tension 



Equation (16 1 is used to compute the surface tension. 



The time t\ is chosen to be when the amplitude a n [given 
by ( 15 1] of the perturbation from a circle exceeds 3 pixels. 
The time £2 is chosen so that for slowly growing patterns 
about 10 images (60 s at 1 frame every 6 s) are collected 
in t<i — ti; for rapidly growing patterns about 20 images 
(10 s at 2 frame/s) are obtained. 

The velocity V was calculated by projecting the lo- 
cal normal to the next later good interface. The use of 
splines allowed for the intersection to occur between in- 
terface points, so that the velocity could be computed 
more precisely. Sums of points of the interface were used 
as approximations of the contour integrals. The wet- 
ting correction (the second integral) contains the capil- 
lary number, Ca = fiV/a, which in our experiments was 
in the range 9 x 10~ 5 <^<2x 10~ 2 . This range of 
Ca is slightly larger than the range that Tabeling and 
Lichaber observed the 2/3 power law relation of the film 
thickness [cf. (pi)] [35]. 



RESULTS 



Harmonic moments 



We have calculated the harmonic moments Mk(t) by 
evaluating the integrals in (16) for 16 bubbles grown in 



a cell with a 125 /im gap and for 10 bubbles grown in a 
Hclc-Shaw cell with a 384 /im gap. For each bubble the 
values of Mk(t) were computed for 50-500 data points 
spaced at time intervals 0.5-6 s. 

We emphasize that our method for determining har- 
monic moments works well for asymmetric as well as 
symmetric bubbles. We chose to attempt to develop sym- 
metric bubbles as they grew in order to track accurately 
the time evolution of particular moments Mk- Hence 
for each bubble the pumping rate was adjusted in real 
time according to (17), as described in Section V.B (cf. 
Fig. |2j. Most of the bubbles evolved toward a targeted 
ri-fold symmetry, which ranged from 5-fold (low pump- 
ing rate) to 14- fold (high pumping rate). However, some 
bubbles had an initial shape that was too irregular to 
evolve into an approximately n-fold bubble during the 
course of the growth; Fig. 4 is an example of such a 
bubble. 

Our main result is that all observed moments Mfc(t) 
(except Mo) decay in time, as Fig. [3] illustrates. The 
next subsection presents results confirming that dM^/dt 
is proportional to the surface tension a (neglecting wet- 
ting correction) [cf. ( 14 1] , in accord with Richardson's 
result that all moments are conserved in the zero surface 
tension case. 




100 200 300 400 

time (s) 

FIG. 3: (a) An air bubble growing in oil in the Hele-Shaw cell 
(plate separation, 125 /im). Adjacent interfaces are separated 
by 56 s, and the maximum span of the bubble is 8 cm. The 
time-dependent pumping rate selects the 7-fold symmetry (cf. 
Fig. [2J. (b) Amplitudes of moments 6, 7, 8, 10 and 14, com- 
puted using note that different moments have different 
units. All moments decay in time, in contrast to the zero 
surface tension case where the moments are constant. 



B. Test of theory and determination of the surface 
tension 



In the theoretical expression for Mkit-i)— M).(ti), equa- 
tion (16 1, all quantities can either be determined from 



our measurements or from independent measurements. 
Hence the theory can be directly tested with no ad- 
justable parameters. 

We chose to present the test of theory as the ratio 
of the surface tension deduced from (16), which we call 
o measured, to the reference value of surface tension mea- 



FIG. 4: An air bubble growing without a clean n-fold sym- 
metry. Adjacent interfaces are separated by 20 s. 



be deduced from (161, where the integrals and velocity 



V are determined from analyses of the bubble patterns, 
while the viscosity ju and gap b are measured indepen- 
dently. 

An example of the results deduced for the surface ten- 
sion ratio is shown in Fig. [5] which was computed for 
the non-symmetric bubble pattern in Fig. [4] For this 
bubble, the surface tension deduced without a wetting 
correction is 30% larger than the reference value at short 
times, where the front velocity is large (see Fig. |4t, while 
at long times the front velocity has become smaller and 
the difference between the ratios with and without the 
wetting correction is only a few percent. Further, both 
ratios are within a few percent of unity, so the value of 
surface tension deduced from theory is equal to the ref- 
erence value within the experimental uncertainty. 

A measure of the symmetry of a bubble is given by 
the spectrum of the harmonic moments, i.e. a plot of the 
normalized moment amplitude, as shown in Fig.[6]for the 
symmetric bubble in Fig. [3] and for the non-symmetric 
bubble in Fig. [4] The approximately 7-fold bubble in 
Fig. [3] has a spectrum with significant components at 
7, 14, and 21 [Fig. 6(a)], while the the non-symmetric 
bubble in Fig. [4] has a spectrum with a larger number 
of significant components, including those at 5, 6, 7, 8 
[Fig-gb)]. 

The mean value for the surface tension a meaS ured was 



deduced from (16 1, including the wetting correction, for 



sured by the Wilhemy Plate method, a 



reference- 



Thus 



surface tension is treated as an unknown whose value can 



each of the 26 bubbles studied, where each bubble was 
evaluated at 50-500 times during the course of its growth; 
this involved the numerical evaluation of the integrals 
in (16 1 for more than 6000 data sets. The result is 
v measured = 18±4 mN/m, where the uncertainty includes 
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FIG. 5: The ratio of the surface tension value deduced from 
theory for the data in Fig. 4 to the value independently mea- 
sured by a traditional method (see text). The upper points 
(x) are without a wetting correction, and the lower points 
(•) are with the wetting correction in ( 16 1 [the term with the 
factor 3.8]. The error bars show the standard deviation for 
each group of points. The value of — ti in ( 16 1 was 10 s. 
The horizontal axis corresponds to t\. 
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FIG. 6: Spectra of the nond imensional moment amplitudes 

J\/ fc |/|M |T~ [equation (15 1] for two bubbles: (a) A bubble 
with approximately 7- fold symmetry [Fig. [3|a)], where the 
main moments are 7, 14, and 21. The moments were calcu- 
lated at the 112th second of the experiment, (b) A bubble 
without n-fold symmetry (Fig. [4b, where there is a broad 
spectrum with no dominant moment. The moments were cal- 
culated at the 64th second of the experiment. 



the surface tension is 23 ± 6 mN/m.] 

The uncertainty in our result for a measure d unfortu- 
nately arises in part from a possible systematic error of 
9%. The uncertainty in gap thickness (about 4%) con- 
tributes significantly to the overall uncertainty. Also, 
an intermittent problem in camera synchronization could 
have introduced timing errors in some of the data. Other 
possible sources of error include the discretization of the 
integrals and the approximations introduced by the the- 
oretical analysis. The range of applicability of the 2/3 
scaling for the film thickness is smaller than the range 
of Ca produced in our experiments. Another intriguing 
possible source of error is the 3.8 factor in the wetting 
correction |34| . which has its origin in a numerical inte- 
gration done in the 1974 PhD Dissertation by Ruschak 
[3"5] . If the numerical factor of 3.8 were changed to 1.5, 
then the mean of the distribution in Fig. [7] would cor- 
respond to unity. Note that the wetting correction be- 
comes small at long times because the growth velocity 
V becomes small if the bubbles are grown according to 
(171. 



The statistics of the results are presented in Fig. [7J as 
histograms of the ratio cr measure d/cr re f erence , calculated 
both with and without the wetting correction. The dis- 
tribution including the wetting correction is narrower be- 
cause the uncorrected data depend on the pattern growth 
velocity. The mean value of the ratio with wetting correc- 
tion is 0.84±0.09 (standard deviation, statistical uncer- 
tainty only); the mean value without wetting correction 
is 1.10±0.1 (standard deviation). 

Future experiments can straightforwardly reduce the 
systematic uncertainty by an order of magnitude to a 
level approaching 1%. Such experiments should focus on 
low order moments at intermediate times, to minimize 
the correction due to wetting and prevent backward mo- 
tion of the interface. 

Bubbles in our experiments were all studied for pump- 
ing rates Q > 0. Alternatively, one could stop the growth 
(set Q=0) and observe the relaxation of a bubble. We 
conducted a few experiments in this way and found that 
the extracted surface tension values were typically 15- 
25% lower than those for growing bubbles. However, 
for relaxing bubbles the correction for the oil wetting 
film, which is re-absorbed as the interface retreats, is not 
known. 



VI. DISCUSSION 



both the statistical uncertainty and the estimated sys- 
tematic uncertainty (see next paragraph). The result 
for Omeasured agrees within the experimental uncertainty 
with <J re f erence — 21.1 ± 0.1 mN/m, determined by the 
Wilhelmy Plate method. Thus the theory is quantita- 
tively confirmed within the experimental uncertainty. [If 
the wetting correction in ( 16 1 is neglected, the result for 



A. Harmonic moments 

We have presented the first demonstration that robust 
harmonic moments Mfc(i) can be deduced from labora- 
tory experiments. The Mk(t) evolve in time because of 
nonzero surface tension [Fig. 3(b)]; otherwise all mo- 
ments except the area M would be conserved. The har- 
monic moments are purely geometric quantities [cf. (|6J> ] 
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FIG. 7: The distribution of the ratios of surface tension de- 
duced from theory, a measured , to er re/erence , determined by 
the Wilhelmy Plate method. The solid line represents com- 
putations that include the wetting term, while the dotted line 
neglects the wetting correction. The bin size is 0.05. A few 
outlying points are outside of the abscissa's range. 



and do not involve approximations, fluid properties, or 
experimental parameters, such as the wetting correction, 
viscosity, or the gap thickness. 

Further, we have obtained an equation for the dynam- 
ics of the harmonic moments, equation (16 1, which in- 



volves quantities that can all be directly determined by 
experiment, thus providing a quantitative test of the the- 
ory of harmonic moments. The result from measurements 
on 26 bubbles for two different Hele-Shaw cell gap thick- 
nesses is that theory and experiment agree within the 
experimental uncertainty (about 20%). This agreement 
was obtained both on nonsymmetric bubbles and on bub- 
bles with approximate n-fold symmetries varying from 5 
to 14. This agreement is robust and, within our experi- 
ments, does not depend on how far a bubble has deviated 
from a circle. Further, we have indicated how it should be 
straightforward to reduce the experimental uncertainty 
by an order of magnitude, to test the theory at about a 
1% level. 

For convenience we have used an expression obtained 
by Bataille (17) to adjust the pumping rate Q(t) to grow 
bubbles that are approximately n-fold symmetric. We 
find that if an n-fold bubble is grown initially from a 
slightly perturbed circle, the n-fold symmetry is retained 
even for large amplitude fingers, far beyond range of ap- 
plicability of (17 1, which was obtained from a linear sta- 
bility analysis. 



B. Implications 

Harmonic moments Mj~ form a complete basis for any 
complicated time domain D, provided that D is analytic 
and singly connected. These moments are exceptionally 
useful for Laplacian growth phenomena such as viscous 
fingering because all other known representations involve 
coefficients that change quickly because of instability, 
thus creating analytic and computational difficulties. In 
contrast, harmonic moments arc free of these problems, 
no matter how much the domain D deviates from a cir- 
cle. Loosely speaking, for harmonic moments the whole 
instability is simply concentrated in the lowest one, Mq 
(the area of D). 

From a fundamental perspective, the moments Mk for 
Laplacian growth are exactly the basis that linearizes and 
"solves" the zero surface tension problem, generating a 
multitude of exact solutions (summarized in [14] ; see also 
references therein), which are impossible to obtain in any 
other basis. In this sense the harmonic moments consti- 
tute the most natural basis for Laplacian growth. Our 
experiments confirm that for nonzero surface tension, the 
derivative of the moments dMk /dt is simply proportional 
to a (neglecting the wetting correction) , so that the zero 
surface tension limit is smoothly approached. This indi- 
cates that the fundamental theoretical results mentioned 
above, which have been obtained for Laplacian growth 
theory for the case of zero surface tension, are relevant 
to real physical systems. 

In conclusion, the implications of the demonstration 
of the harmonic moments description of viscous fingering 
extends into various branches of mathematical physics 
(cf. Section II) because of the deep connection between 
Laplacian growth and other problems. 
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